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We extend a recently discovered, non-singular 6 dimensional brane, solution to -D = 4 + n di- 
mensions. As with the previous 6D solution the present solution provides a gravitational trapping 
mechanism for fields of spin 0, i, 1 and 2. There is an important distinction between 2 extra dimen- 
sions and n extra dimensions that makes this more than a trivial extension. In contrast to gravity 
in n > 2 dimensions, gravity in n = 2 dimensions is conformally flat. The stress-energy tensor 
required by this solution has reasonable physically properties, and for n = 2 and n = 3 can be made 
to asymptotically go to zero as one moves away from the brane. 



I. INTRODUCTION 



The brane world models introduced in 0, Q have provided a new framework for addressing various questions in 
particle physics, such as the hierarchy problem, using non-compact extra dimensions (see 131,11 IE Umll for earlier 
^sj ' work on non-compact extra dimensions and on solutions in extra dimensions). In these models one must give a 
mechanism for trapping or localizing fields with various spins (0, ^, 1, 2) on the 4 dimensional brane of the observable 
' Universe. The condition for determining whether a multi-dimensional field is localized or not is if the integral of its 
Lagrangian over the extra coordinates is finite. In the original 5 dimensional models it was found that it was not 



. possible to gravitational trap all types of fields with the same warp factor. Spin and spin 2 fields were trapped on 

■ the brane with a decreasing, exponential, warp factor, while spin i fields were trapped with an increasing warp factor 
^f-^ , 9] . Spin 1 fields were not localized at all |0| . In the 6 dimensional case with two extra dimensions slightly better 

trapping behavior was found 0, 1 and 2 fields were localized with a decreasing, exponential warp factor but spin 
, i fields were again localized with an increasing, exponential warp factor. In Il2|| it was shown that one could localize 

■ the zero modes of all these fields with an increasing, non-exponential warp factor in 6 dimensions with one extra 
time-like dimension. In 0, it was shown that all these fields could be localized in 6 dimensions with both extra 
dimensions being space-like. In particular [l^ flSi] gave an exact, non-singular, non-exponential trapping solution. 
Given the improving localization behavior with increasing spacetime dimension the natural inclination would be to 

D look for trapping solutions with still more extra dimensions. However, it was pointed out in 16] that there is an 
important distinction in going from n = 2 to n > 2 extra dimensions. Gravity in 2 dimensions is conformally flat. 
^ , For spacetime with greater than 2 dimensions vacuum spacetime can have a non-zero Weyl tensor. The present work 
shows that it is possible to extend (with some restrictions) the non-singular, trapping solution of [Hill toD = 4 



n 



[ dimensions. (Earlier work on the localization of gravity in n > 2 extra dimensions with higher curvature terms can 
^ > be found in 0|). The solution localizes all types of spin fields with a non-exponential, warp factor. The draw back 
■ ■ ■ ' of the solution is that the stress-energy tensor is selected or tuned so as to give these good features of the solutions. 

However the stress-energy tensor is physically reasonable in that its magnitude is peaked on the brane and decreases 
as one moves away from the brane. It will be shown that for the n — 2 and n = 3 cases it is possible to have the 
stress-energy tensor go to zero as one moves away from the brane. For rt > 4 one can have part, but not all of the 
stress-energy tensor asymptotically go to zero. 

It is possible to localize all spin fields in the 5D and 6D models with exponential warp factors if one introduces 
non-gravitational interactions. However, for reasons of simplicity, and since gravity couples universally to all types of 
matter, it is preferable to find a purely gravitational trapping mechanism. The solution presented here does give such 
a pure gravitational trapping mechanism for all types of spin fields. The warp factor is non-exponential in contrast 
to the work in pi .2. ■!&■ .IQj 
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II. D-DIMENSIONAL EINSTEIN EQUATIONS 



The system we consider is D — 4 + n dimensional gravity with a cosmological constant A, and some matter fields. 
The action is 



]\fn+2 

^— (^i? + 2A) + L„ 



(1) 



where \/ —^g is the determinant, M is the fundamental mass scale, ^ R is the scalar curvature, A is the cosmological 
constant and L is the Lagrangian of matter fields. All of these quantities are D = 4+n dimensional. The D-dimcnsional 
Einstein equations resulting from this action arc 

""Rab - \gAB ^i? - A.g^B + J^Tab , (2) 

Capital Latin indices run over A, B, ... = 0, 1, 2, 3, D. The stress-energy tensor, Tab, is determined from the matter 
Lagrangian, Lm, in the standard way. 

The ansatz we will use is a simple extension of refs. [Tsl IT^ with two scale functions - one for the 4-dimensional 
brane and one for the n extra spatial dimensions. 

ds^ = (j,^{r)r]a,fi{x'')dx"dx'^ - X{r){dr^ + r^dni_^) . (3) 

The metric of ordinary 4-dimensional spacetime, rjapix'^), has the signature — , — , — ). The Greek indices a, f3, ... ~ 
0,1,2,3 refer to the coordinates of these 4-dimensions. It is assumed that ansatz Q only depends on the extra 
coordinate r through the 4-dimensional conformal factors, 0^, and A. df2^_]^ is the solid angle for the (n — 1) sphere 
and is explicitly given by 

dnl_^ = del_^ + sin^ er,-idel_^ + sin^ 6i„_i sin^ er,-2del_j, + ■ • ■ + n ^i"^ ^''^^i ('^) 

1=2 

The angles have the following ranges: Q < 6i < 2tt and < 6*^ < tt for i = 2;...,n — 1. The radial coordinate has the 
range < r < oo. 

The above generalization of the 6D ansatz of [13. IT^ has an important difference from the generalization of the 6D 
metric ansatz considered in [TH ITslIT^ . In these latter papers the scale factor, A(r), multiples only the angular part 
of the metric 



= 4?{r)r]ap{x'')dx°'dxl^ - dr"^ - \{r)de^ , (5) 

whereas in 0,^3 the scale factor multiplies the entire extra dimensional part. Both the 6D version of ^ and Q 
are locally the same. This can be seen by defining dr of (O via X(r)dr'^ = df^ where f is the radial coordinate of 
(O . However there is a global distinction between Q and . The metric in jSJ corresponds to a cone-like geometry 
such as occurs with string-like defect. For example, if in |SJl one has A(r) = (1 — n)^r^ then the geometry of the extra 
space will exhibit an angular deficit of 2?™. A similar A(r) in Q would not result in an angular deficit. 

For the 6D ansatz of |13L and its generalization to 4-|-n given in Q the metric functions, A(r) and <^(r), are 
conformal factors for an n-dimensional Euclidean metric and 4-dimensional Minkowski metric respectively. Thus as 
long as the metric functions are well behaved (i.e. not equal to or oo) the metric will be non-singular. Since we are 
looking for non-singular solutions, we choose the metric to take the form 

The nonzero components of the stress-energy tensor Tab are taken to be of the form 

Tp^i, = -g^^F{r), = ~gijK{r), T^, = . (6) 

The two source functions, F{r) and K{r), depend only on the radial coordinate r. 

As a final simplification we assume that 4-dimensional metric is fiat {riap{x'^) — rjap) and that the 4-dimensional 
cosmological constant is zero. This yields 



^-R/ii/ — -r]fj,u*R — , 



(7) 
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where ^Rap, '^R are the 4-dimensional Ricci tensor and scalar constructed from r^n^ With these simplifications and 
the ansatze of ^ and lO the D-dimensional Einstein field equations ^ become |23| 



+ in - 1) 



r2A 



r2A 



12- 



+ (n-l) 



(/) r^A 



r2A 



4 2 



(/) (/) A 



12- 



+ (n-2) 



'■'(r^A)' (r^A)" n-5/(r2A) 



r^A 



r2A 



(n-2) 



lA'(r2A)' n-2 
" 2TT2A ~ 

lV (r2A)' n-3 
2 A r2A r2 

= 2A(^-A 

M"+2 



(8) 
(9) 



(10) 



where the prime = 3/ dr. These equations are for the aa, rr, and 99 components respectively. 

Subtracting the rr component of (|5J| from the 99 component of 1)11) |l to cancel the source and cosmological terms 
yields 



\ cf r^J + M A 2\\) rX 



. 



(11) 



For n — 2 one has the simplification that the second term above does not appear. Following j2flj we look for solutions 
in which each term in parentheses above is zero separately 



(f) \ <f> r(j) 



A solution to the first equation in (|12|l is 



A(r) 



A 2 I A / r\ 



(12) 



(13) 



where p is an integration constant with dimensions of length. In [20| it was shown that the trapping solution of 



could be generalized to more than two extra dimensions. The following metric functions 



(r) 



A(r) 



2 „2 



4ac p- 
(r2 + c?Y 



r'',c'' 



(14) 



We will find 



solve H12f) . with a and c constants. In the solution was more general in that one had r'^,c 

a similar selection of the exponent when we generalize the non-singular solution of [l^ll5j| to n > 2 dimensions. The 
solution 1)14(1 becomes singular at r = c and thus was taken as a solution for r > e where e is the brane width, and 
the condition, e > c, holds. Outside r — e the scale factor for 4-dimensional space, (^^(r), asymptotically approaches 
the value a^. If the brane width, e, is small then the 4 dimensional scale factor changes very quickly from 0^(0) = 1 
to <^2^oo) — a^. In 6D the solution in 1)14)) is supported by a stress-energy tensor whose source ansatz functions, 
F{r), K{r) where proportional to The solution (j){r) in 1)14)) was shown ]l3.l20j | to localize fields of all spins up to 
2 on the brane. 

In mill a solution was given which had similar properties to the solution of 1)14)) (non-exponential warp factor 
which localized fields of all spins up to spin 2), but which was non-singular. Here it is shown that, with some 
modifications and restrictions, it is possible to extend this solution to n extra spatial dimensions. This is more than 
a trivial exercise since for n = 2 the extra space is conformally flat, while for n > 2 this is not necessarily the case. 
Also for both n — 2 and n = 3 it is possible to make a choice of parameters {i.e. e, a. A) such that F{r) and K{r) 
both have the nice feature of going to zero as r — > 00. For rt > 4 one can make either F[r) or K{r) go to zero, but 
not both. 

It is straightforward to check that the non-singular solution of [l3l [Tsj l 



c + ar^ 
(p- -\- r2 



A(r) 



2(a - I)c2p2 
(c2 -I- r2)2 



(15) 
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satisfies both equations from (|12|l : a and c are constants. The solution given in was more general than (|15|l having 
the form 



(16) 



where the exponent of r could take values other than 2 and was restricted by |6| > 2. In the present case if one inserts 
the 4'{r) from into the second equation in H12|l (the first is automatically satisfied via the choice in one 
obtains 



3 /A'" ' 



2 V A / r\ 2r2 



(17) 



which only vanishes for the choice h = ±2. Thus for n > 2 the b = ±2 solution of |0| is selected. In [3| the b > 2 
solutions gave a scale factor, </'(r), which monotonically increased as one moved away from the brane. The b < 2 
solutions gave a scale factor, 0(r), which decreased monotonically as one moved away from the brane. For the solution 
in IjlSI) the increasing (f>{r) corresponds to a > 1 and the decreasing (f>{r) to a < 1. It was pointed out in jl^ that 
both increasing and decreasing solutions give localization of fields with spins ranging from to 2. 

Some of the parameters in (|15|1 can be fixed. The width of the brane, e, is given by the inflection point of (f){r), i.e. 
(j)"{r = e) = 0. This conditions gives c — 3e^. Next we require that A(r = 0) = 1, i.e. that on the brane one has a 
6-dimensional Minkowski metric. This gives — 2(^a-i) ■ Putting these parameters back into (I15II yields 

3e2 + ar2 

3e2 + r-2 ' ^'■^ " (3e2 + r2)2 

Now taking (p and A from the above equation and inserting them into equations JHJ and ^ we obtain expressions for 
the source functions F{r) and K{r) 



F{r) _^ ^ -6e2n(n-3a + 2) ^ 2 (3(n + 2) + 3a2(n + 2) - 2a(n2 + 2n + 6)) r2 



M"+2 (3e2 + ar2)2 (3^2 +^^2)2 

-2an (-3 + a{n + 2)) r'^ 
^ 3e2(3e2 + ar2)2 

K{r) _^ ^ -6e2(n-l)(n-4a + 2) ^ 4 (2(n + 2) + 2a2(n + 2) - a(n2 + n + 10)) r2 



M"+2 (3e2_^^^2)2 (3£2_|_„^2)2 

-2a(n - 1) (-4 + a{n + 2)) 
^ 3e2(3e2 + ar2)2 



(19) 



(20) 



For n = 2 and using the relationship = "^^"^ from TS*! one can show that F{r), K{r) reduce to the expressions in 
It is straightforward to check that the last equation (|10l) is also satisfied by (/)(r), A(r) from (|18|l and K{r) from 
(|20|l . One can also check that the stress-energy tensor defined by F{r) and if(r) above satisfy energy-momentum 
conservation given by 

V^Tab = -^=dA{V^9T^'') + rgzji^'''' = (21) 
V- 5 

In terms of the source ansatz functions this yields 

K' + A—{K-F)^0 . (22) 

Using the F{r) and K{r) from ifT!?)) one finds that indeed is satisfied. Thus taken together ((TH|l ((T^ (I^OJ) give a 
non-singular, trapping solution to the 4 -t- n dimensional Einstein equations © © (|10|l . 



III. LOCALIZATION OF ZERO MODES 



As with the 6D solution of [TsLITsI l the present solution provides a gravitational trapping for fields of various spins: 
0, ^, 1, 2. The condition for a field to be trapped on the brane is that the integral of the Lagrangian over the extra 
coordinates is finite. 



(23) 
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where Ls is the matter Lagrangian for fields with spin s = 0, ^, 1, 2. We will outline the trapping of the zero modes 
for each of these types of fields by applying the solution (|18|l to the analysis given in |23|. The Lagrangians for the 
various spin fields will contribute different factors of 0, A and r to the common factor of y'-'^g = (/)4A"/2r.«-i in 
For a scalar field the specific from of the action is 

So = -lJ d'^x^/^g g^'^dM^dN^ (24) 

The equations of motion from this scalar Lagrangian have a zero-mode solution 0| of the form <^{x'^^) = 4>{x^)uq, 
with Mo being the constant solution for the extra dimensions, and 0(x'^) satisfies the Klein-Gordon equation on the 
brane rj^^'^ d^du(t>{x) = 0. Thus the scalar action becomes 

So = -f|^«o rfr02^"/2r"-i J d^x r;^"^ 9^09,0 (25) 

The standard gamma function above comes from the integral over dfln-i namely J dfl„-i — ■ The extra factor 
of comes from g'^'^ in the Lagrangian. Inserting the solution into the radial integral we obtain 



2"+i I r(2±li) ^ ' ' r{2 + n) 



which is finite. 

The action for the spinor case is 



Si/2 = / d'^xy/^ ^iT^'Dm-^. (27) 



1/2 

The curved spacetime gamma matrices (F*^) are related to the flat spacetime gamma matrices (7*^) via the inverse 
vielbein ejg, F*'^ = e^7^. The vielbein is defined via gMN = ^m^n'Hmn- Dm'^ — {dhi + j^^^f Inp)'^ is the covariant 
derivative. The spin connection oj^/' is defined in terms of the veilbein and its derivatives in the standard way. In 
[20| the zero mode for the 4 + n dimensional Dirac equation was found to be of the form '^{x^) = ^{x^^)u{r)x{Oi), 
where ip{x'^) satisfied the massless 4 dimensional Dirac equation, 'y'^df^ip = 0, and xi^i) satisfied the angular part of 
the higher dimensional Dirac equation. The specific forms of -0 and x were not needed in the analysis of the trapping 
of the fields. The u{r) part of the higher dimensional spinor was found 2^ to have the form 

_ 71-1 

uir) = Ci/2r' (A^r) ' (28) 

where C1/2 is an integration constant. In contrast to the other fields, for the spinor field the extra space part of the 
wavefunction is not simply a constant. Inserting all this into the spinor action in H27|) yields 

Si/2 = f dr^b^'X^r^'-'u^r) ^ d^!„_lx'(^^^) J d^xibtj''d^^ 

drc^-^X^ ^dnn^lxHO^) I d^xi>ird^^ (29) 



The extra factor of in the first exp ression comes from the veilbein of the higher dimensional gamma matrix, F'^. 
The integral over dQ,n-i is finite |2^ and the integral over r is 

r drr'X- = r dr^^^^ - ^ (30) 
Jo Jo 3e^ + ar^ 2^ 

which is also finite. Thus the integral over the extra coordinates is finite and the fermions are localized. 
For the vector, spin 1 case the action is given by 

Si^-\l d^xV^g g''''g'''FMRFMs , (3i) 
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where Fmn = QmAm — O^Am- There is a zero mode solution [23| of the form Am{x'^) = a^(a;'')ito with uq the 
constant solution over the n extra space dimensions and satisfies the normal 4 dimensional vacuum Maxwell 
equation, d^{d^ay — d^afj) — d^f^„ — 0. Inserting all this along with the solution of (fH^ into (pTT|l gives 



Si = -^^4 I dr\-/\-' I d^xrj^'^v'y.xf,^ (32) 



The radial integral is evaluated as 



oo 



n^2n„n-l qn/2 n rz ~r f n\ 

Since the radial integral is finite the vector fields are also localized on the brane by this solution. 

The analysis for the spin-2 is similar to the spin-0 case and results in the same radial integral (see ^3)- Thus 
the present solution also localizes the 4-dimensional spin-2 graviton on the brane, since the integral 1)26(1 is finite. 
As for the 6D solution, the present solution localizes spin-0, ^, f , 2 onto the brane via gravity alone. In all cases the 
integral over the extra coordinates is finite independent of whether the 4D scale function, 0^, is increasing (a > 1) or 
decreasing (a < 1). Thus, as pointed out in both types of solutions localize the fields. This is connected with 
the fact that the extra space scale function, A, is independent of a, and decreases to zero away from the brane in all 
cases. 



IV. CONCLUSIONS 



The above results are in some sense expected since the scale factors, 0(r), A(r), given in H18I) arc of the same form 
as in the 6D case in However in contrast to the general 6D solution found in 15], the present higher dimensional 
solution selects out the h = ±2 exponent for the dependence of </>(?'). 

The behavior of the ansatz functions, F{r) and K(r), depends on the various parameters, n, e, A and a. In particular 
A and a play the major role in determining the values of F{r) and K{r) at r = and r = 00. For a large choice of 
parameters which lead to 5-like behavior for F{r) and K{r) (i.e. peaked near r = and decreasing as r ^ 00). Let 
us now focus on the r oo limit of F{r) and K{r). It is possible to have either F{r) or K{r) or both go to zero for 
some choices of parameters, a, A, e, and n. To see this we first write the asymptotic values of F{7') and K{r) from 

Fir ^00) = AM"^^+^"^^""^(^"/^" + ^)) 
Kir 00) = AM"- + 2in-l)M-^^iA ain + 2)) 



It is possible to make i^(oo) = by choosing 

6n 



An + 2n2 - 3e2A 

and it is possible to make if (00) = by choosing 

8(n-l) 



2(n-l)(n + 2)-3e2A 



(35) 



(36) 



If A = in H35|l and H36|) then a = and a — respectively. Thus for A = one can only have Fir) or Kir) 
asymptotically go to zero with a < 1 solutions. One must have A > if one wants Fir) or Kir) asymptotically to go 
to zero and have a > 1. It is also possible to have both Fir) and Kir) go to zero at r ^ cx). Setting H35|l equal to 
((36|l and solving for ie^K gives 



2. 2n(7i + 2)(n-l) 



3,-A = ^ '-^ '- (37) 

n — 4 

Inserting this into either H35() or H36() gives a = |^ which yields a < 1 for n > 2. For the 6D case considered in 
the focus was on a > 1 solutions and thus there it was not possible to have both F(r) and Kir) asymptotically go 
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to zero. Thus there are two possibihties to have both F{qo) = and K{oo) = 0: n = 2 and n = 3 which correspond 
to a = i and a — ^. For n > 4 and a < and the solution in (|15|l becomes singular. The n = 4 case also has 
problems. First the scale function 0(r) becomes singular at r = cx). Also from (|37|l one finds that either e or A or 
both are infinite. Finally, if one recomputes the source ansatz functions, F{r),K(r) in this case one finds that they 
always go to oo at r = oo. However an interesting result is that for n = 3 one has a non-singular branc solution with 
sources that approach zero away from the brane. In 16] it was argued that for local defects (i.e. solutions where 
the stress-energy was strictly zero outside the brane or decreased exponentially) it might not be possible to construct 
non-singular solutions for n > 2. Our n = 3 solution evades this restrictions since the stress-energy goes to zero as a 
power law rather than exponentially or being strictly zero outside the brane. 

In this paper we have extended the non-singular, 6 dimensional brane, world trapping solution of [l^ Il5l | to 
n > 2 extra dimensions. This is more than a simple mathematical exercise since, as pointed out in pla |. there is 
a significant difference between n = 2 and n > 2 extra dimensions. As with the 6D solution, the present solution 
provides a universal, gravitational trapping mechanism for fields of various spins from spin to spin 2. Considering 
the magnitude of the stress-energy ansatz functions necessary to realize this trapping solution one finds physically 
reasonable properties for a range of parameters n, A, e and a. In addition for n = 2 and n = 3 it is possible to have 
the entire stress-energy tensor go to zero at r = oo. For n > 4 one can have either F{r) or K{r) go to zero, but not 
both. We have not addressed the question of the stability of this brane solution in the present paper. 
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